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Scattering by regular black holes: Planar massless scalar waves impinging upon a 

Bardeen black hole 
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Singularities are common features of general relativity black holes. However, within general 
relativity, one can construct black holes that present no singularities. These regular black hole 
solutions can be achieved by, for instance, relaxing one of the energy conditions on the stress-energy 
tensor sourcing the black hole. Some regular black hole solutions were found in the context of 
nonlinear electrodynamics, the Bardeen black hole being the first one proposed. In this paper, we 
consider a planar massless scalar wave scattered by a Bardeen black hole. We compare the scattering 
cross section computed using a partial-wave description with the classical geodesic scattering of a 
stream of null geodesics, as well as with the semiclassical glory approximation. We obtain that, for 
some values of the corresponding black hole charge, the scattering cross section of a Bardeen black 
hole has a similar interference pattern to a Reissner-Nordstrom black hole. 

PACS numbers: 04.70.-s, 04.70.Bw, 11.80.-m, 04.30.Nk 


I. INTRODUCTION 

Black holes (BHs) are among the most interesting ob¬ 
jects of general relativity (GR). Although GR is a highly 
nonlinear theory, BHs come out of it with a very simple 
structure. Standard GR BH solutions are parameter¬ 
ized by their mass, charge and angular momentum [1] 
(see, e.g., Refs. [2, 3] for interesting counterexamples of 
the previous statement). Although related to the earliest 
predictions of GR, the strong field regime of BHs is still 
an experimental challenge [4-6]. Notwithstanding, the 
observational data presently available suggest that BHs 
populate basically all the galaxies in the Universe [7]. 

Although very successful in explaining the available 
data, standard GR BHs suffer from one of the main prob¬ 
lems of GR: the presence of singularities. Hawking and 
Penrose indeed showed that, for some hypotheses on the 
gravitational collapse, the formation of singularities in 
BHs would be unavoidable [8, 9]. These singularities 
were conjectured by Penrose to be hidden by a horizon 
[10, 11], and were claimed to be possibly avoided within 
an improved theory of gravity (extension or modification 
of GR) [12]. 

The study of BHs without singularities can help us 
to understand the role played by singularities in astro¬ 
physics. Still within GR, one can obtain BHs with¬ 
out singularities — dubbed regular BHs — by relaxing 
one of the energy conditions on the stress-energy ten¬ 
sor. Bardeen proposed the first regular BH solution [13], 
which was later identified as a solution for a nonlinear 
magnetic monopole [14]. Since then, other regular BHs 
appeared in the literature, in different scenarios (see, e.g., 
Refs. [15-17], and references therein). Moreover, regular 
BHs can be relevant in the context of quantum gravity 
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[18, 19], and some of them reproduce the quantum weak 
field regime of GR [20]. 

Studies of scattering by BHs have been extensively 
made [21], Reference [22] presented the results of the 
scattering of all basic massless (spin 0, 1/2, 1 and 2) 
fields by Schwarzschild BHs. Moreover, the shadows of 
BHs [23, 24] may become visible with future telescopes, 
like the Event Horizon Telescope [25] , and the scattering 
of light, considering wave and semiclassical approxima¬ 
tions, may be important in anticipating subtle charac¬ 
teristics of the shadows. However, a careful study of the 
scattering of fields by regular BHs is still lacking in the 
literature [26]. 

The line element of the Bardeen BH can be written as 
ds 2 = -f{r)dt 2 + f(r)~ * 1 dr 2 + r 2 (d0 2 + sin 2 9d<p 2 ), (1) 
where the lapse function f(r) is given by 

2 Mr 2 

Hr) = 1 - ( 2 ) 

The Bardeen BH has a structure similar to the Reissner- 
Nordstrom (RN) BH (see, e.g., Ref. [15]). For q < q ext = 
4M/(3-\/3) (henceforth, without loss of generality, we 
shall assume q > 0) the spacetime has two horizons and 
for q = q ext the horizons degenerate, characterizing the 
extremal case. Following Ref. [27], we shall display our 
results in terms of the normalized charge Q = q/q ex t [28]. 

It is interesting to note that the variation of the charge 
of the regular BH changes considerably the deflection an¬ 
gle of light rays [29] . In Fig. 1 we plot null geodesics com¬ 
ing from infinity with an impact parameter of b = 5.2 M, 
for different values of the BH charge. The behavior is 
qualitatively similar to the case of RN BHs [30]. We see 
that, adjusting the BH charge, we can have scattering in 
basically any direction. 

The remainder of this paper is organized as follows: 
In Sec. II we review the main aspects of the classical 
geodesic scattering and semiclassical glory approximation 
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that we also consider the case in which the null geodesic 
rotates (one or many times) around the BH before going 
to infinity (for instance, see the solid curve of Fig. 1). It 
is interesting to note that the classical scattering formula 
given by Eq. (3) describes very well the planar-wave case 
for small scattering angles, although it gives discrepant 
results for moderate-to-high scattering angles, as we shall 
see in Sec. IV. 

Let us now obtain b(x) through a geodesic analysis. 
Without loss of generality, we shall restrict the geodesic 
motion to the plane 9 = 7t/2. From the line element (1), 
we can write, for null geodesics 
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f{l/u}u 2 , 


( 4 ) 


FIG. 1. Geodesics approaching a Bardeen BH from infinity 
with an impact parameter of b = 5.2 M, for different values 
of the BH charge. The Schwarzschild case (Q = 0) is also 
exhibited (solid line). 


where we have defined u = b = L/E is the impact 
parameter given in terms of the constants of motion 

E = ft and L = r 2 ip, (5) 


to compute the differential scattering cross section. In 
Sec. Ill we present the partial-wave method used to com¬ 
pute the scattering cross section of planar massless scalar 
waves. In Sec. IV we present the results for the scat¬ 
tering of planar massless scalar waves impinging upon a 
Bardeen regular BH, comparing the three different ap¬ 
proaches used to compute the scattering cross section. 
In Sec. V we end up with our final remarks. Throughout 
this work we use G = c = H = 1 and metric signature 

( — > +> +j +)• 

II. CLASSICAL SCATTERING AND 
SEMICLASSICAL GLORY 

In this section we investigate the scattering by BHs by 
using two approaches: classical geodesic scattering and 
the semiclassical glory approximation. These approaches 
allow us to foresee some of the aspects of the scattering 
cross section obtained within the full partial-wave anal¬ 
ysis. 


A. Geodesic scattering 


and the overdot denotes differentiation with respect to 
an affine parameter of the curve. 

Differentiating Eq. (4) with respect to (p, we obtain 


d 2 u 

dp 2 


u 2 df{l/u) 
2 du 


- uf(l/u). 


( 6 ) 


Solving Eq. (6) with the appropriate boundary condi¬ 
tions, one can obtain the geodesics followed by mass¬ 
less particles, such as the ones shown in Fig. 1. The 
smallest positive root of the right-hand side of Eq. (6) 
corresponds to the radius of the critical orbit for null 
geodesics, u c = l/r c . Substituting its value in the right- 
hand side of Eq. (4) and setting it to zero, we obtain 
the impact parameter associated with the critical orbits, 
b c . Going in the other way around, i.e., choosing a value 
b > 6 C , the smallest root of the right-hand side of Eq. (4) 
is Uq = 1/ro, where r 0 is the turning point — or the 
radius of maximum approximation — of the geodesic. 

Finally, by integrating Eq. (4) in the case of scattered 
geodesics, we obtain 


u 0 



o 


f(l/u)u 2 


- 1/2 


du. 


( 7 ) 


The analysis of null geodesics in the Bardeen spacetime 
can be seen in Ref. [27]. For the classical approximation 
of the scattering we may consider a stream of parallel null 
geodesics coming from infinity. In this case, the analysis 
of Ref. [31] suits the problem of classical scattering by 
Bardeen BHs. The classical scattering cross section is 
given by 


da 

dtt 


sm x 




db(x) 

dx 


( 3 ) 


where 6(x) is the impact parameter associated with a 
scattering angle %. The summation in Eq. (3) is such 


The deflection angle following directly from Eq. (7) is 
given by 


0(6) = 2a(6) - 7r. (8) 

By inverting Eq. (8), one obtains 6(0), and substituting 
it in Eq. (3), one obtains the classical scattering cross 
section. Plots of the classical scattering cross section ob¬ 
tained by Eq. (3) are exhibited in Fig. 2, where we 
compare the Bardeen cases with the Schwarzschild one. 
We see from Fig. 2 that the increase of the BH charge 
contributes to an increase of the BH classical scattering 
cross section. 
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FIG. 2. Left panel: Classical scattering cross section of Bardeen BHs, with Q = 0.5, 0.8 and 1, and for the Schwarzschild BH 
(Q = 0). Right panel: Classical scattering cross section of Bardeen BHs, normalized by the Schwarzschild case. 


B. Glory scattering 

The interference that occurs between scattered partial 
waves with different angular momenta is not taken into 
account by the classical formula (3). In order to obtain 
a scattering cross section that takes into account the in¬ 
terference processes, we need to perform a wave analysis 
of the problem. Before going into the full wave analy¬ 
sis, however, it is interesting to apply an approximate 
method that works remarkably well for high scattering 
angles ( 9 ~ 7r) and that captures some key features of the 
scattering cross section in this regime, including the in¬ 
terference process: the semiclassical glory approximation 
[32]. Indeed, one of the main advantages of this semi¬ 
classical approximation is that one can find an analytical 
formula that gives some physical insight for the width 
of interference fringes in the scattering cross section as 
well as the intensity of the scattered flux for 9 ~ n. As 
a semiclassical approximation, it is valid for will 1, 
although it can still reproduce remarkably well some re¬ 
sults for ojM ~ 1. 

The semiclassical formula for the glory scattering by 
spherically symmetric BHs is given by [32] 

J| s (w6 s sin 0), (9) 

6 — 7T 

where b g is the impact parameter of backscattered rays 
(9 = 7r), J 2 s{x) is the Bessel function of the first kind (of 
order 2s), and s is the wave spin. In our case, since we 
are considering a scalar wave, s = 0. We note that there 
are multiple values of b g corresponding to the multiple 
values of the deflection angle, namely 0 = n + 2nn , with 
n = 0,1, 2, that result on back-scattered rays. All 
the rays scattered close to the the backward direction 
( 9 ~ 7r) contribute to the glory scattering, but the most 
important contribution comes from the n = 0 case. The 
next contribution, n = 1, has an intensity that is about 
0.2 % of the n = 0 one in the Schwarzschild case, and 


about 0.8 % in the case of the extreme Bardeen BH. This 
is a consequence of the derivative \db/d9\g =n in Eq. (9) 
getting rapidly suppressed as n increases. In fact, the 
values of b g for rays that pass multiple times around the 
BH are very close to each other and also to b c . Here, 
we consider only the most important contribution to the 
glory scattering. 

Once we have the knowledge of the glory scatter¬ 
ing formula, Eq. (9), we only need to determine b g 
and \db/d9\g =n in order to obtain the glory scattering 
cross section. Therefore, we apply Newton’s method 
and numerical integration to obtain the parameters b g 
and \db/d9\g =n . Numerical results for r c , b c , b g and 
bg\db/d9\g= n are presented in Fig. 3. From these re¬ 
sults and Eq. (9), we may expect that (i) interference 
fringes get wider and (ii) backscattered flux intensity is 
enhanced for higher values of the BH charge. Expecta¬ 
tion (i) comes from the fact that the interference fringe 
width is inversely proportional to b g , as indicated by the 
argument of the Bessel function in Eq. (9). Moreover, 
Fig. 3 shows that b g decreases monotonically as Q in¬ 
creases. We also note from Eq. (9) that the scattering 
intensity is proportional to b 2 g \db / d9\g =n , and to the wave 
frequency. As shown in Fig. 3, although b g decreases 
monotonically with the increase of Q , b g \db/d9\g =n in¬ 
creases monotonically with Q, justifying expectation (ii). 

The above analysis may be compared with the results 
for the glory scattering from RN BHs [30]. In the latter 
case, b g decreases monotonically with the increase of the 
BH charge, while \db/d9\g= n increases monotonically (c/. 
Fig. 9 of Ref. [30]). Therefore, considering the change 
of the BH charge, the behavior of the parameters b g and 
\db/d9\g— n are qualitatively the same for Bardeen and RN 
BHs. In the case of RN BHs, however, the combination 
bg\db/d9\g =v does not increase monotonically with the 
charge — as it happens for Bardeen BHs. Instead, the 
glory scattering amplitude as a function of Q presents a 
local minimum in the RN case (cf. Fig. 10 of Ref. [30]). 




db 

d9 
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FIG. 3. Glory scattering parameters for Bardeen BHs with 
varying charge, considering only the dominant contribution 
(0 = 7r). We note that, with the exception of \db/dO\$= n , 
all the important parameters related to the glory scattering 
decrease monotonically with the increase of Q. 


In Sec. IV, we compare results obtained from Eq. (9) 
with partial-wave results, exhibiting excellent agreement 
in the regime 8 < n. 


III. PLANAR WAVE SCATTERING 

Planar massless scalar waves, represented by the wave 
function <E>, are described by the Klein-Gordon equation 

~^=!d a (V^99 ab d b $) = 0 . ( 10 ) 

Here we shall be interested in monochromatic plane 
waves. We have 

^ = E (11) 

Im 

where Y l m {8, ip) are the scalar spherical harmonics. Sub¬ 
stituting Eq. (11) into Eq. (10), we obtain the following 
radial equation: 


+ V 4>( r ) ~ w2 ) ^( r ) = °’ ( 12 ) 

in which r* is the tortoise coordinate, defined through 
dr* = /(r) _1 dr, and 


V^,(r) = / 


Kl + 1) , 

r 2 ' 



(13) 


is the scalar field effective potential. The scalar field 
potential is localized, going to zero at both asymptotic 
limits of r* (infinity and horizon) [27]. 

Plane waves coming from infinity can be described in 
terms of the so-called in modes. These modes are purely 


incoming from the past null infinity, obeying the follow¬ 
ing boundary conditions: 


(b(r) ~ / ^ 1 d" ^uZ-Rj) as r* —> +oo (r —> +oo), 

^ ~ \TuiRii, as r* -+ -oo (r -+ r h ), 

with 

A A U) 

R I = e ~ iur -Y l -jr> ( 15 ) 

3=0 T 
N 

i?// = e-^*^(r-r ft )M«, (16) 

3=0 

where \TZuji \ 2 and \Tu>i \ 2 are the reflection and transmis¬ 
sion coefficients, respectively. Flux conservation implies 
that \Ti u i \ 2 + \Tuji\ 2 = 1. Note that the summations in 
Eqs. (15) and (16) are required to keep track of the 
convergence of the solutions. The numerical infinity and 
horizon are chosen such that V^(r) <C w 2 at the bound¬ 
aries. 

The scalar differential scattering cross section for 
Bardeen BHs can be written in terms of partial waves 
as [21] 

I - W")' 2 . < 17 > 

where 


+ !) [ e2i5i(w) ~ Pilose) (18) 


1=0 


is the scattering amplitude, with the phase shifts {Si) 
given by 

e 2iSl (") = {-l) l+l n u i. (19) 


IV. RESULTS 

In order to obtain the phase shifts to compute the 
scattering cross section via the partial-wave method, we 
have applied a fourth-fifth Runge-Kutta method to solve 
the radial equation (12). We have typically started with 
the near-horizon condition at r s = l.OOOlr/,, and the 
outer boundary (numerical infinity) chosen depends on 
the value of l. Results were obtained with boundary 
conditions (14), as well as with alternative conditions in 
terms of spherical Hankel functions (see, e.g., Eq. (18) 
of Ref. [30]). Both conditions lead basically to the same 
results. Since the sum in Eq. (18) does not converge 
very quickly, because of the Coulomb characteristic of 
the problem, we have applied the convergence method 
first introduced by Yennie et al. [33] , and first applied to 
the BH scattering problem by Dolan et al. in Ref. [34]. 

In Fig. 4 we show the scattering cross section for 
Bardeen BHs with different charges (Q = 0.5, 0.8,1), as 
well as for the Schwarzschild BH, and tfw = 2. We 
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FIG. 4. Scattering cross sections for Bardeen BHs considering 
different BH charges. We also plot the Schwarzschild case, for 
comparison. We see that the value of the charge affects the 
fringe widths, while the (avarage) amount of scattered flux 
remains basically the same. 

see that the fringe widths increase with the increase of 
the BH charge. This, as anticipated by the semiclassi- 
cal analysis of Sec. IIB, is in accordance with the fact 
that b g decreases monotonically with the increase of Q , 
as previously seen in Fig. 3. The amount of scattered flux 
(on average) remains basically the same. These general 
behaviors are similar to the ones presented by the RN 
BHs [30], 

Figure 5 presents comparisons of the numerical scalar 
scattering cross sections for Bardeen BHs with the ap¬ 
proximated geodesic and glory results. We see that the 
glory results fit remarkably well the numerical results for 
large angles (9 < 7r), while the geodesic results fit well 
the small-angle region. This very good agreement can 
also be considered as a consistency check of our results. 

The glory approximation can be used to capture most 
of the features of the back-scattered wave. Some caution, 
however, should be taken when one considers the glory 
intensity in the backward direction (normalized by the 
wave frequency). To illustrate this, in Fig. 6 we plot the 
amplitudes of the back-scattered wave, for Mlo = 1, 2 
and 3, computed through the partial-wave method and 
through the glory approximation, as a function of the 
Bardeen BH charge. We see that the results obtained 
via the partial-wave method oscillate around the one ob¬ 
tained using the glory approximation. This agrees with 
the analysis presented in Ref. [30] for RN BHs. 

In Fig. 7 we compare the differential scattering cross 
sections of Bardeen and RN BHs. While RN and Bardeen 
BHs with the correspondent charge produce different 
scattering patterns - illustrated by the top-left panel 
of Fig. 7 — we can have configurations with different 
charges that produce almost the same scattering pattern. 
A similarity of the patterns also happens in the case of 
absorption cross sections, when the critical impact pa¬ 
rameter (b c ) of the RN and Bardeen cases are the same 


[27]. Here, however, the similarity of the scattering cross 
sections intensifies when b g for the RN and Bardeen cases 
match. The similarities are illustrated in the top-right 
and bottom panels of Fig. 7, where we show the scatter¬ 
ing cross section for a RN BH with Q = 0.753 and for a 
Bardeen BH with Q = 1, for different values of the fre¬ 
quency. The scattering flux intensities are different for 
intermediate-to-high scattering angles, while the inter¬ 
ference widths are essentially the same for all scattering 
angles. 


V. FINAL REMARKS 

We have computed the scalar scattering cross section of 
regular Bardeen BHs. Numerical results were compared 
with both geodesic and glory approximations, and we 
have found excellent agreement within the validity limits 
of each approximation. 

From the glory approximation, it is known that the 
interference fringe widths depend inversely on the im¬ 
pact parameter of backscattered waves, b g . The classical 
analysis from geodesics shows that b g decreases monoton¬ 
ically with the increase of Q. Therefore, we expect that 
the interference fringes get wider as Q increases. This 
was confirmed from our numerical results obtained via 
partial-wave method. 

Comparison of Bardeen BHs with RN BHs reveals that 
the scattering of these two kinds of BHs can be similar 
but not identical. By similar we mean that in both cases 
the scattering cross section presents (i) intense oscilla¬ 
tions in the near-backward scattering, (ii) rapidly grow¬ 
ing flux amount and smoother oscillations for smaller an¬ 
gles, and (iii) similar results for very small scattering an¬ 
gles. (i) is a consequence of the strong interference be¬ 
tween rays passing by the opposite sides of the BH, as 
it is well described by the glory approximation in both 
cases, (ii) is a consequence of the fact that for small scat¬ 
tering angles both b and \db/dd\ increase as 6 diminishes, 
and the difference between paths followed by neighboring 
rays becomes smaller, weakening interference effects. We 
can conclude that the main contribution to the scatter¬ 
ing cross section for very small angles comes from rays 
with high impact parameters [35]. We may treat such 
cases in the weak-field regime, where the main contribu¬ 
tion to the gravitational interaction comes from the BH 
mass, i.e., f[r) ~ 1 — 2 M/r + 0{r~ n ), where n = 2 in the 
case of RN BHs, and n = 3 in the case of Bardeen BHs. 
This explains (iii), i.e., why, for BHs with the same mass, 
in the regime of small angles, all results tend to be the 
same, regardless of the nature and value of their charge. 

The results presented in this paper reinforce in a way 
the results presented in Ref. [27], implying the conclusion 
that some properties of Bardeen BHs can be very similar 
to those of RN BHs (with different charge). In this sense, 
we conclude that it may be difficult to discriminate regu¬ 
lar BHs from the standard ones, as far as absorption and 
scattering of scalar plane waves are concerned. It should 
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FIG. 5. Comparison of partial-wave, semiclassical glory and classical geodesic approaches for the differential scattering cross 
section, for Muj = 3 (in the first two cases), and different values of the Bardeen BH charge [Q = 0.5 (left) and Q = 1 (right)]. 
The semiclassical glory approximation reproduces very well the results for backscattered waves (9 ~ tt), while the classical 
approach works well for small scattering angles. 



Q 


FIG. 6. Glory intensity in the backward direction (normalized 
by the wave frequency) (9 = 7r) as a function of the Bardeen 
BH charge. We see that the amplitude computed through 
the partial-wave method oscillates around the one computed 
through the glory approximation. 


be interesting to extend the analyses presented here and 
in Ref. [27] to the scattering and absorption of waves with 
higher spins and compare with recently obtained results 
for RN BHs [36-40], 
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